Superstripes and the excitation spectrum of a spin-orbit-coupled Bose-Einstein 

condensate 
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Using Bogoliubov theory we calculate the excitation spectrum of a spinor Bose-Einstein condensed 
gas with equal Rashba and Dresselhaus spin-orbit coupling in the stripe phase. The emergence 
of a double gapless band is pointed out as a key signature of Bose-Einstein condensation and 
of the spontaneous breaking of translational invariance symmetry. In the long wavelength limit 
the lower and upper branches exhibit, respectively, a clear spin and density nature. For wave 
vectors close to the first Brillouin zone the lower branch acquires an important density character 
responsible for the divergent behavior of the structure factor and of the static response function. 
The sound velocities are calculated as functions of the Raman coupling for excitations propagating 
orthogonal and parallel to the stripes. Our predictions provide new perspectives for the identification 
of supersolid phenomena in ultracold atomic gases. 
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The search for supcrsolidity represents a challenging 
topic of research in different areas of condensed matter 
and atomic physics (for a recent review see, for exam- 
ple, [l|). Supersolidity was first predicted in the pio- 
neering works by Andreev and Lifschitz Q, Leggett Q, 
and Chester It is characterized by the co-existence 
of two spontaneously broken symmetries. The breaking 
of gauge symmetry gives rise to off-diagonal long range 
order yielding superfluidity, while the breaking of trans- 
lational invariance yields diagonal long range order char- 
acterizing the crystalline structure. First experimental 
efforts towards the search of supersolidity were carried 
out in solid helium The strongly interacting nature 
of this system makes however the effects due to Bose- 
Einstein condensation extremely small and no conclusive 
proof of supersolidity is still available in such a system 
|6|. More recently, systematic attempts to predict the 
occurrence of a supersolid phase have been carried out in 
atomic gases with dipolar [7H9( and soft core, finite range 
interactions However, these configurations have 

not yet been experimentally realized in the quantum de- 
generate phase required to observe the new effects. 

The recent realization of spinor Bose-Einstein con- 
densed gases with spin-orbit coupling 17-2fjj is open- 
ing new perspectives in the field. In systems with equal 
Rashba and Dresselhaus coupling and for small values 
of the Raman coupling theory in fact predicts the oc- 
currence of a stripe phase where translational invariance 
is spontaneously broken [2l| - |23j . Actually these systems 
are periodic only in one direction and can be considered 
as superfluid nematic liquid crystals. Experiments are al- 
ready available in the relevant range of parameters, but 
no direct evidence of the density modulations is still avail- 
able, due to the smallness of the contrast and the micro- 
scopic distance separating consecutive stripes. A phase 
transition has been nevertheless detected [20( at values 
of the Raman coupling below which theory predicts the 



occurrence of the stripe phase. 

The purpose of the present work is to show that the 
excitation spectrum of the gas in the stripe phase ex- 
hibits typical supersolid features, like the occurrence of 
two gapless bands and the divergent behavior of the static 
structure factor for wave vectors approaching the bound- 
ary of the Brillouin zone. The excitation spectrum is 
measurable in Bragg spectroscopy experiments, so the 
experimental characterization of the new phase should 
not represent a major difficulty. 

Spin-orbit-coupled Bose-Einstein condensed gases can 
be described using the mean-field Gross-Pitaevskii pic- 
ture. The interaction is zero ranged and is characterized 
by the values of the scattering lengths associated with 
the two hyperfine states involved in the Raman process 
(we will limit here the discussion to spinor Bose gases). 
This differs from the case of other systems, like dipolar 
gases, where the origin of the supersolid phase is associ- 
ated with the finite range of the force [IH . The validity of 
the Gross-Pitaevskii approach can be tested a posteriori 
by evaluating the quantum depletion of the condensate. 

We will use the following single-particle Hamiltonian 
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accounting for the effect of two counter propagating and 
polarized laser fields, where fco is fixed by the momentum 
transfer of the two lasers, while is the Raman coupling, 
accounting for the intensity of the laser beams causing 
the transition between the two spin states. The occur- 
rence of the term 5a z /2 has been ignored in h , since 
we will consider situation where the effective magnetic 
field 8 is zero (experimentally this can be achieved with a 
proper detuning of the two laser fields). The Hamiltonian 
([I]) can be formally derived by applying a unitary trans- 
formation to the Hamiltonian in the laboratory frame de- 
scribing the system in the presence of two detuned, spin 
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polarized laser fields [26[. The unitary transformation 
consists of a local rotation in spin space around the z- 
axis, which causes the appearance of the spin-orbit term 
proportional to p x <J z . 

Remarkable properties of this Hamiltonian are its 
translational invariance and, for values of the Raman 
coupling fl smaller than 2fcg, the occurrence of a dou- 
ble minimum structure in the single-particle energy at 
momenta p x = ±ki with ki = fccy! — ^ 2 /4^0' capable 
to host a Bose-Einstein condensate. This structure is at 
the origin of new intriguing features, like the existence 
of a spin polarized plane- wave phase and of an unpolar- 
ized stripe phase | 



22j at even smaller values of fl, 



resulting from the spontaneous breaking of translational 
symmetry. 

The emergence of the stripe phase, which is the main 
object of the present Letter, is the consequence of the 
competition between the density and spin-density inter- 
action terms in the mean-field Hamiltonian 

H int = \Jd 3 r [(g + g n ) n(vf + (g - g n ) S (r) 2 ] (2) 

where n(r) = n.f(r) + nj.(r) and s(r) = nf(r) — nj.(r) 
correspond to the total and spin densities, while g a p = 
4ira a /3 (a, ft = t, 4) are the coupling constants in the 
different spin channels, a a p being the corresponding s- 
wave scattering lengths. For simplicity, we have assumed 
that the intraspecies scattering lengths are equal: = 
a-il [13] ■ The stripe phase emerges only for < g, a 
condition yielding an unpolarized uniform ground state 
in the absence of the Raman coupling f2. It is associated 
with the macroscopic occupation of a single-particle state 
given by the combinations of spinor states of the form 
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where k\ — 7r / ' d is related to the period d of the stripes, 
K = 2nfci, with n = 0, ±1, • • • , are the reciprocal lat- 
tice vectors while a_ kl+ x and b_ kl+ x are expansion co- 
efficients to be determined, together with the value of 
k\ , through a procedure of energy minimization, includ- 
ing the single-particle ([I]) and the interaction ([2]) terms 
in the Hamiltonian. In the stripe phase energy mini- 
mization gives rise to the presence of terms with op- 
posite phase (e ±lklX , e ±3lfcia: , •••) characterized by the 
symmetry condition a_ fel+i j = 6* causing the van- 
ishing of the average spin polarization. The occurrence 
of terms with different phases is at the origin of the den- 
sity modulations. In the limit of weak interactions the 
spinor wave function ([3]) is characterized by the values 
a±k 1 = b T ki = \/ (kg ± fci)/2fco predicted by the single- 
particle Hamiltonian (fTJ), while the higher order expan- 
sion coefficients become vanishingly small. The other 
phases predicted by theory at larger values of fl are the 
plane-wave and the zero momentum phases. In these 



phases only the term e lklX (or e~ iklX ) is present in the 
sum ((3|) with ki ^ in the plane-wave phase and k± = 
in the zero momentum one. 

As already pointed out, in the stripe phase the spin po- 
larization is vanishing but the density modulations cause 
an energy cost in the density term of Eq. ([2]) with respect 
to a uniform configuration. Viceversa in the plane-wave 
phase the density is uniform but the spin polarization 
differs from zero and is given by ki/ka [22j. The final 
energetic balance between the two phases is determined 
by the value of the Raman coupling. In the limit of small 
coupling constants, defined by the condition G\, Gi <C 
k$, where G t = n (g + g n ) / 4 , G 2 = ft {g - gn) /4 and 
n is the average density, the critical value for the Ra- 
man frequency is given by fi cr = 2*^27/(1 + 27) with 
7 = G2/G1 independent of the density, and the stripe 
phase is ensured for values of fl smaller than Q cr . In 



the available experiments with 87 Rb atoms 2^, 28| the 
value of G2 (and hence f2 cr ) is very small. In order to en- 
large the range of values of fl compatible with the stripe 
phase it is then necessary to choose G2 as large as pos- 
sible. This is crucial to produce a significant contrast 
in the density profile which is proportional to Q/Aiq [22T |. 
In the following we will use the values Gi/fc 2 = 0.3 and 
Gijkl = 0.08 yielding O cr /fcg ~ 1.3. In Fig. Q] we show 
the ground state density profile calculated at fi/fcg = 1. 

In order to evaluate the elementary excitations we ap- 
ply Bogoliubov theory by writing the deviations of the 
order parameter with respect to equilibrium as 



-ifit 



u t (r) 



vl(r) 



(4) 

and solving the corresponding linearized time-dependent 
Gross-Pitaevskii equations. The equations are conve- 
niently solved by expanding ^^(r) and «-[-. j,(r) in the 




FIG. 1: Density profile along the x direction. The parameters 
are fi/fcg = 1.0, Gi/ftg = 0.3, and G 2 /fcg = 0.08. 
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FIG. 2: (Color online) Lowest four excitation bands along the 0.5 1 1 .5 2 

x direction. The parameters are the same as in Fig. [T] ^1 



Bloch form in terms of the reciprocal lattice vectors: 



« qta (r) =e- ik ^Y,U^ i& e i **+ iR * (5) 

K 

w qtt ^(r) = ( --^l„. (6) 

K 

where q is the wave vector of the excitation. The 
same ansatz for ipf t ± can be also used to calculate the 
density and spin-density dynamic response function, by 
adding to the Hamiltonian a perturbation proportional 
to e t(i-r-"t)+ n t and ^(q-r-^H^ resp ectively. 

The excitation spectrum predicted by the Hamiltonian 
(HJ)-© has been already calculated in the plan e- wave 
phase as well as in the zero momentum phase |26f where, 
despite the spinor nature of the system, only one gap- 
less branch is predicted as a consequence of the presence 
of the Raman coupling tt. A peculiar feature exhibited 
by the plane-wave phase is also the emergence of a ro- 
tonic structure whose gap becomes smaller and smaller 
as one decreases the value of fi, providing the onset of 
the transition to the stripe phase. 

The results for the dispersion law of the elementary 
excitations in the stripe phase are reported in Fig. [2] for 
the same parameters used in Fig. Q] We have consid- 
ered excitations propagating in the x direction orthog- 
onal to the stripes and labeled with the wave vector 
q x . A peculiar feature, which distinguishes the stripe 
phase from the other uniform phases, is the occurrence 
of two gaplcss bands. At small wave vectors we find 
that the lower branch is basically a spin excitation, 
while the upper branch is a density mode, as clearly re- 
vealed by Fig. [3] (a) where we show the contributions of 
the two gapless branches to the static structure factor 

S(q x ) = N- 1 Y,e\(°\Pq*\ty\ 2 where £ is the band index 
and p qx = J2i e tqxXi is the density operator. A two- 
photon Bragg scattering experiment with laser frequen- 
cies far from resonance, being sensitive to the density 
response, will consequently excite only the upper branch 



FIG. 3: (Color online) Density (a) and spin-density (b) static 
structure factor as a function of q x (blue solid line). The 
contributions of the first (red dashed line) and second (black 
dash-dotted line) bands are also shown. The parameters are 
the same as in Fig. [T] 



at small q x . Bragg scattering experiments actually mea- 
sure the imaginary part of the response function, a quan- 
tity which, at enough low temperature, can be identified 
with the T = value of the dynamic structure factor 
S(fe,w) = N- 1 Ee\(Q\PqJt)\ 2 S(T^ ~ ^ojwhere E m is 
the excitation energy of the i-t\x state [30j|. The spin 
nature of the lower branch is clearly revealed by Fig. [3] 
(b) where we report the contributions arising from the 
two gapless branches to the spin static structure factor 
S a (q x ) = N- 1 ^ |(0|sj^| 2 where s qx = J2 z a zl e^ x > is 
the spin-density operator. Notice that, differently from 
S(q x ), the total spin structure factor does not vanish as 
q x — > 0, being affected by the higher energy bands as a 
consequence of the Raman term in the spin-orbit Hamil- 
tonian (UJ) . The lower branch actually exhibits a hybrid 
nature and when one approaches the Brillouin wave vec- 
tor qB = 2fci it is responsible for the divergent behavior 
of the density static structure factor [see Fig. [3] (a)], a 
typical feature exhibited by crystals. Its direct obser- 
vation would represent an unambiguous evidence of the 
stripe phase. 

The divergent behavior near the Brillouin zone is even 
more pronounced (see Fig. if one investigates the 
static response function x(Qx) = 2 J2e I (0|Pg* K) \ 2 / E eo, 
proportional to the inverse energy weighted moment of 
the dynamic structure factor. The divergent behaviors 
of S(q x ) and ~x(q x ) can be rigorously proven using the 
Bogoliubov [29( and the Uncertainty Principle 3l| in- 
equalities applied to systems with a spontaneously bro- 
ken symmetry. These inequalities are based, respec- 
tively, on the relationships m-i(F)mi{G) > \([F, G])\ 2 
and m (F)m (G) > \([F, G])\ 2 involving the p-th mo- 
ments m p (0) - (|(0|d^)| 2 + |(O|0t|^)| 2 ) ( Ee - E )p 
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FIG. 4: (Color online) Static response as a function of q x (blue 
solid line). The contributions of the first (red dashed line) and 
second (black dash-dotted line) bands are also shown. The 
parameters are the same as in Fig. [T] 



of the l-th strengths of the operators F = J^j e iq:rXj and 
G = T,j(Pxj e- i ^"- qB ^ + H.c.)/2. The commutator 
( [F, G] ) = q x N(e lqsX ), entering the right hand side of 
the inequalities, coincides with the relevant crystalline 
order parameter and is proportional to the density mod- 
ulations of the stripes. The moments m-i(F) and mo{F) 
are instead proportional to the static response x{Qx) and 
to the static structure factor S(q x ), respectively. It is not 
difficult to show that the moments m\{G) and mo(G) are 
proportional, respectively, to (q x — qs) 2 and to \q x — qs\ 
as q x — > qs as a consequence of the translational invari- 
ance of the Hamiltonian, thereby causing the divergent 
behaviors S(q x ) oc l/\q x - q B \ and xfe) oc l/(q x - q B ) 2 
exhibited, respectively, by the static structure factor and 
by the static response. 

In addition to the excitations propagating along x (lon- 
gitudinal modes) another class of Bogoliubov modes is 
predicted in the transverse directions, parallel to the 
stripes. These modes are excited by the density and 
spin-density operators ^2 i e lqvVi and ^2 i cr Z ie lqyyi . Re- 
markably, also in the transverse channel Bogoliubov the- 
ory predicts the occurrence of two gapless spectra. The 
lowest branch has a spin character at small q y . Viceversa 
the upper branch is essentially a density wave, similarly 
to what happens in the longitudinal channel. 

In Fig. [5] we compare the sound velocities of the two 
gapless branches in the longitudinal and transverse di- 
rections. In the transverse direction the density sound 
velocity is well reproduced by the Bogoliubov expression 
^/2Gl (equal to 0.78 k in our case) for the sound veloc- 
ity. The sound velocity in the longitudinal direction is 
always smaller than in the transverse case, reflecting the 
inertia of the flow caused by the presence of the stripes. 
Notice also that the spin sound velocity becomes lower 
and lower as the Raman frequency increases, approach- 
ing the transition to the plane- wave phase. Notice that 
the Bogoliubov solutions in the stripe phase exist also 
for values of f2 larger than the critical value Q cr = 1.3kg, 
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FIG. 5: (Color online) Sound velocities in the first (red) and 
second (black) bands along the x (solid lines) and perpendicu- 
lar (dashed lines) direction as a function of £7. The blue dash- 
dotted line represents the transition from the stripe phase to 
the plane- wave phase. The value of the parameters Gi/fcj) 
and G2 / k$ are the same as in Fig. [1] 



due to the first order nature of the transition. 

We have finally checked that the quantum depletion of 
the condensate, due to the fluctuations associated with 
the Bogoliubov solutions, are always very small, thereby 
confirming the validity of the mean field approach. 

In conclusion we have shown that the excitation spec- 
trum in the stripe phase of a spin-orbit-coupled Bose- 
Einstcin condensate exhibits a double gapless band struc- 
ture, typical of supersolids. We predict that at small 
wave vectors the lower and upper branches have, respec- 
tively, a spin and density nature. The lower branch, 
whose gapless nature is the consequence of the breaking 
of translational symmetry, is responsible for the divergent 
behavior of the static structure factor as the wave vector 
approaches the border of the Brillouin zone. The ex- 
perimental verification of the new dynamic features pre- 
dicted in this Letter is expected to provide a significant 
advance in our understanding of systems exhibiting si- 
multaneously off-diagonal and diagonal long range order. 
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